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We demonstrate that manganite heterostructures of the form 

{Insulator) / {LaMnOi) n / {SrMnOz)n / {Insulator) show strong multiferroicity with magnetic 
ordering being induced by ferroelectric charge ordering. Ferroelectric polarization is produced by 
charges leaking from LaMnO-j, side to the SrMnO-j, side on account of the effective nearest-neighbor 
repulsion of electrons resulting from electron-phonon interaction. The origin of the ferromagnetism 
is due to double exchange when a particle virtually hops to its unoccupied neighboring site and 
back. Furthermore, we also demonstrate that a striking giant magnetoelectric effect occurs due to 
the strong coupling between the charge and spin orderings. 
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I. INTRODUCTION 

Oxides such as manganites display a rich interplay of 
various orbital, charge, and spin orders when rare-earth 
dopants are added to the parent oxide. While significant 
progress has been made in characterizing the bulk doped 
materials, the heterostructures produced from two par- 
ent oxides is an area of active research!. Several novel 
phenomena, that have no counter part in the bulk sam- 
ples, seem to emerge from the intcrface/heterostructure 
physics 2 -. In fact, the challenge is to technologically ex- 
ploit this new physics and develop new useful devices that 
can meet future challenges such as miniaturization, dis- 
sipationless operation/manipulation (such as read/write 
capability), energy storage, etc. 

There has been a revival in multifcrroic research partly 
due to improved technology, discovery of new compounds 
(such as YMnOj,, T6M712O5), need for devices with 
strong magnetoelectric effect, eto 3 -. Majority of the mul- 
tiferroics studied are bulk materials where it is not yet 
clear why the magnetic and electric polarizations coexist 
poorly^. With the advent of improved molecular-beam- 
cpitaxy technology one can now grow oxide heterostruc- 
tures and explore the possibility of strong multifcrroic 
phenomena as well as large interplay between ferroelec- 
tricity and magnetic polarization. There are numerous 
mechanisms for magnetoelectric effect reviews for which 
can be found in Refs. la-0- At the interface of a magnetic 
oxide and a ferroelectric/dielectric oxide, magnetoelectric 
effect of electronic origin has been predicted by some re- 
searchers. Upon application of an external electric field, 
not only the magnitude of moments can be changed^, 
but in some cases the very nature of magnetic ordering 
can be changed^. 

In this paper we predict a giant magnetoelectric effect, 
not at the interface, but away from it, in an oxide-oxide 
heterostructure (see Fig. 1). Our objective is to present 



a plausible/possible unnoticed multiferroic phenomenon 
in manganite heterostructures and point out the asso- 
ciated large/striking magnetoelectric effect. Jahn- Teller 
electron-phonon interaction is shown to be key to under- 
standing both multiferroicity and magnetoelectric effect 
in our oxide heterostructure. Here, we exploit the fact 
that manganites have various competing phases that are 
close in energy and that by using an external pertur- 
bation (such as an electric or a magnetic field) the sys- 
tem can be induced to alter its phase. We show that 
there is a charge redistribution (with a net electric dipolc 
moment perpendicular to the interface) and a concomi- 
tant ferromagnetism due to the optimization produced 
by the following two competing effects: (i) energy cost 
to produce holes on the LaMnOz (LMO) side and ex- 
cess electrons on the SrMnOz (SMO) side; and (ii) en- 
ergy gain due to electron-hole attraction (or electron- 
electron repulsion) on nearest-neighbor Mn-sites induced 
by electron-phonon interaction. The charge polarization 
is akin to that of a pn-junction in semiconductors al- 
though the governing equations arc different. The key 
results of our analysis arc: (i) the interface charge den- 
sity at the LMO-SMO interface is 0.5 electrons/site and 
the charge in the layers perpendicular to the interface 
rapidly attain the bulk value as one moves away from 
the interface due to the fact that the dielectric constant 
is not large; (ii) the LMO-SMO interface is ferromagnetic 
and persists to be ferromagnetic for a few layers adjacent 
to the interface on both sides of the interface; (iii) fur- 
thermore, since fcrroelectricity and ferromagnetism have 
a common origin [i.e., holes (electrons) on LMO (SMO) 
side], there is a striking interplay between these two po- 
larizations leading to a giant magnetoelectric effect. 

The rest of the paper is organized as follows. In Sec. II 
we introduce our phenomenological Hamiltonian (based 
on cooperative Jahn- Teller physics) using which we de- 
duce the charge distribution in Sec. Ill and establish 



the ferroelectric effect of our oxide-oxide heterostructure. 
Next, in Sec. IV we show how the charge profile across 
the interface produces magnetic phases (much like the 
phase diagram of LSMO) and show how the interface 
charge affects magnetism on both sides of it. We close in 
Sec. V with our concluding observations. 



II. POLARONIC HAMILTONIAN 

In our (Insulator) / (LaMn0 3 ) n / (SrMn0 3 ) n / (Insulator) 
heterostructure depicted in Fig. 1, we expect to realize 
the entire phase diagram of La\^ x Sr x Mn0 3 starting 
with the x = phase at the Insulator-LMO interface 
and sequentially evolving to the x = 1 phase at the other 
end involving SMO-Insulator interface. Since, majority 
of the LSMO phase diagram (including the end regions 
near x = and x = 1) is taken up by insulating phases, 
we expect that there is no effective transport in the 
direction normal to the oxide-oxide interface. Further- 
more, we will show that the calculated concentrations 
in any layer on both the LMO and SMO sides turn out 
to be such that the bulk Lai- x Sr x Mn0 3 (LSMO) at 
this concentration is always in the insulating regime. 
Thus, we consider a system of localized electrons: This 
is also in tune with the insulating behavior reported 
in (SrMn0 3 ) n / (LaMn0 3 ) 2n super lattices!^. Addi- 
tionally, at strong electron-phonon coupling, one can 
see that band width is significantly diminished. Other 
reasons for localization, besides narrow band electrons 
being localized by disorder (such as interface roughness) , 
could be anti-ferromagnetic order on the SMO side and 
cooperative lattice Jahn- Teller distortion on the LMO 
side. Then, for analyzing the charge distribution normal 
to the interface, the starting Hamiltonian is assumed to 
have the following phcnomenological form: 
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FIG. 1. (Color online) Schematic showing the 

(Insulator) / (LaMnOi)N / (SrMnOi)N / (Insulator-) het- 
erostructure. Each of the labelled N layers on both LMO 
and SMO sides, as well as in the Interface, comprise of 
manganese-oxide (MO) layers. 
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where the first coefficient j^ g 2 cjQ is due to electron- 
phonon interaction and represents nearest-neighbor 
electron-electron repulsion brought about by incompat- 
ible distortions of nearest-neighbor oxygen cages sur- 
rounding occupied Mn ions. The pre-factor 7g p can de- 
pend on the phase - for instance, in the C-type antifer- 
romagnetic (AFM) regime of LSMO, jl p is expected to 
be large because occupancy of neighboring d z i orbitals 
is inhibited in the z-direction; while in the A-type AFM 
regime (corresponding to undoped LaMn0 3 ), 7* is ex- 
pected to be weaker because of compatible Jahn- Teller 
distortions on neighboring sites. Here we will assume for 
simplicity that 7* is concentration independent and that 
0-1 < 7c P < 1- Next, the coefficient j 2 p t 2 /(g 2 uj ) results 
from processes involving hopping to nearest-neighbor and 
back and is present even when we consider the simpler 
Holstcin modelii or the Hubbard-Holstcin models. The 
pre-factor 7^ varies between 1/2 (for non-cooperative 



electron-phonon interaction) and 1/4 (since for coopera- 
tive breathing mode in one-dimensional chains j 2 =1/3, 
which should be more than in C-chains)^. Now, since 
manganites are two-band systems, t 2 /g 2 uo is an effec- 
tive term which has contributions from hopping between 
various orbitals on neighboring sites and the polaronic 
energy associated with different orbital occupancy. Fur- 
thermore, although nj is the total number in both or- 
bitals on site j, it can only take a maximum value of 
1 due to strong on-site electron-electron repulsion and 
strong Hund's coupling. Next, to make the above Hamil- 
tonian furthermore relevant for manganites, one needs 
to consider Hund's coupling between core spins and itin- 
erant electrons. This leads to invoking the double ex- 
change mechanism for transport. Then, the hopping 
term iy between sites i and j in Eq. ([1]) is modified to 
i(l-n j+S ), (1) be t ijy /Q.5[l + (Si -Sj/S 2 )]. Lastly, we ignored the ef- 
fect of superexchange interaction in the starting effective 
Hamiltonian in Eq. ([T]) because its energy scale is signif- 
icantly smaller that the polaronic energy term g 2 ujQ and 
thus its effect on charge distribution can be neglected. 



III. CHARGE PROFILE 

For a localized system, we only need to minimize the 
interaction energy which is a functional of the electronic 
density profile. The Coulombic interaction energy result- 
ing from the electrons leaking from the LMO side to the 
SMO side, is taken into account by ascribing an effec- 
tive charge +1 (hole) to the LMO unit cell (centered at 
the Mn site) that has donated an e g electron from the 
Mn site and effective charge —1 (electron) to the SMO 
unit cell (centered at the Mn site) that has accepted 
an e g electron at the Mn site. The Coulombic energy 
that results is due to interactions between these ±1 ef- 



fective charges. The net positive charge on the LMO side 
and net negative charge on the SMO side will produce a 
charge polarization (or inversion asymmetry). Since the 
ferroelectric dipole is expected to be in the direction per- 
pendicular to the oxide-oxide (LMO-SMO) interface, we 
assume that the density is uniform in each layer for cal- 
culating the density profile as a function of distance z 
from the insulator-oxide interface. The Coulombic inter- 
action energy per unit area due to leaked charges is the 
same for both LMO and SMO regions and is given, in 
the continuum approximation, to be 



E m , = ±f dzD(zf 



(2) 



where D(z) is the electric displacement and is given by 



D(z) = ± 



-eE R 



dy4irep(y) 



(3) 



with + (-) sign for LMO (SMO) side. Furthermore, 
p(z) is the charge density, e the charge of a hole, e the 
dielectric constant, E ea:t an external electric field along 
the z-direction, and L the thickness of the LMO (SMO) 
layers. 

The ground state energy per unit area [corresponding 
to the effective Hamiltonian of Eq. (JTJ] can be written, 
in the continuum approximation, as a functional of p(z) 
for both LMO and SMO as follows: 
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(4) 

where C = 6 is the coordination number and a is the lat- 
tice constant. In arriving at the above equation we have 
approximated rii+s + rii-s ~ 2nj. Furthermore, for the 
ground state, we expect i«j\/0.5[l + (Si ■ S 3 /S 2 )] = tij, 
i.e., minority charges completely polarize the neighboring 
majority charges. 

We will now minimize the total energy given below 

Erotai = 2E po i + 2E cou i, (5) 

by setting the functional derivative SETotai/Sp(z) =0. 
This leads to the following equation 







-C 1 [l-2a 3 p(z)]+2C 2 / dy 
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27re 2 /e; and 



E ext = eE ea; t/(47re). The above equation, upon taking 
double derivative with respect to z, yields 



da 



* d 2 p{z) 

dz 2 



C 2 p(z) = 0. 



(7) 



The above second-order differential Eq. (J7J and Eq. © 
admit the solution 



Pi?) 



1 cosh(mz) 
2a 3 cosh(mL) 



rnE e 



sinh[m(L 



cosh(mL) 



(8) 



where m = -J C 2 / (Cia 3 ) . It is important to note that 
for the manganese-oxide (MO) layer at the LMO-SMO 
interface (i.e., at z = L), the density is 0.5 electrons/site 
and that it is independent of the applied external electric 
field! Furthermore, for realistic values of the parameters, 
the charge density rapidly changes as we move away 
from the oxide-oxide interface (i.e., after only a few 
layers from the interface) and attains values close to the 
bulk value. Lastly, as required for zero values of field 
Eextj w e get p(0) — s- when L — > oo. Thus, although 
we used the continuum approximation, our obtained 
density profile is qualitatively realistic as it has the 
desired values at the extremes z = L and z = with 
the density away from the LMO-SMO interface rapidly 
falling for not too large values of e. 



For our calculations we used the following values for 
the parameters: a — 4 A; e = 20; 0.05 cV < ujq < 
0.07 eV; 2 < g < 3; and 0.2 eV < t < 0.6 cV. The density 
profiles for both LMO and SMO sides depend only on 
e, C\ and Fi ext . We observe from Figs. 2 and 3 that 
for smaller values of C\, as distance from the interface 
increases, the density falls more rapidly while the density 
change due to electric field rises faster. 




LMO/SMO 



FIG. 2. (Color online) Minority-charge density p 

in various layers (on both LMO and SMO sides) 
of a (Insulator) / (LaMnO^)^ / (SrMnOs)^ / (Insulator) het- 
erostructure for a — 4 A, e = 20, and C\ — 0.19. Dis- 
played results are for (a) p at E e xt — kV/cm and (b) 
Ap = p(E ext = 100 kV/cm) - p(E ext = kV/cm). This 
system does not show a significant magnetoelectric effect due 
to low density p away from the interface. 




crostructurc is expected to be zero except in the vicinity 
of the interface where (due to proximity effect) it will 
be FM and can be modeled using a percolating picture. 
Given the above scenario, as can be expected, we find 
that the FM region on the LMO side can be drastically 
enhanced (at the expense of the A-type AFM region) by 
an electric field inducing holes on the LMO side. On the 
hand, the electric field has only a small effect on the per- 
colating FM cluster that is adjacent to the oxide-oxide 
interface and is on the SMO side. 



A. SMO side 

We will first consider the SMO side and show that the 
magnetization decays as we move away from the LMO- 
SMO interface. We derive the largest FM domain as fol- 
lows. On account of nearest-neighbor repulsion (as given 
in Eq. (TT])), the interface (which is half-filled) has e g elec- 
trons on alternate sites. In fact, to minimize the intcr- 



FIG. 3. (Color online) Minority-charge density p 

and magnetization M/fj,B in various manganese- 
oxide layers (on both LMO and SMO sides) of a 
(Insulator) /(LaMn03) 3 /(SrMn03)'i/ (Insulator) het- 

erostructure for a — 4 A, e = 20, and C\ = 0.93. 
Depicted figures are for (a) p(E ext — kV/cm); (b) 
Ap = p(E ext = 100 kV/cm) - p(E ext = kV/cm); and 
(c) M/p B at E ext = kV/cm and E ex t = 100 kV/cm. 
Manganese-oxide (MO) layer 1 on the LMO side undergoes 
spin reversal when an external field E ex t is applied. 



IV. MAGNETIZATION DISTRIBUTION 

We will now obtain the magnetization for a het- 
erostructure by considering its lattice structure unlike 
the case for the density profile where a continuum ap- 
proximation was made. Thus we can take into account 
the possibility of antifcrromagnctic order besides being 
able to consider ferromagnetic order. 

First, we observe that the bulk Lai- x Sr x MnOs (be- 
low the magnetic transition temperatures) is A-type an- 
tifcrromagnctic (AFM) for < x < 0.1 ; ferromagnetic 
insulator (FMI) for 0.1 < x < 0.15; and ferromagnetic 
(FM) metal for 0.15 < x < 0.5i£. Hence, we model 
the LMO side of the heterostructure as a ferromagnet in 
regions with hole concentrations that are large but less 
than 0.5 with the magnetic structure abruptly switch- 
ing to an A-type AFM as the hole concentration values 
decrease to about 0.1 or less. Next, we note that for 
0-5 $ x < 1.0, the LSMO bulk system is always an 
AFM insulator with A-type AFM for 0.5 < x < 0.65; 
C-type AFM for 0.65 < x < 0.95; and G-typc AFM 
for 0.95 < a; < 1.0. Thus, from a magnetism point-of- 
view, the magnetic moment on the SMO side of the het- 
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FIG. 4. (Color online) Minority-charge density p and 
magnetization M/ps in various manganese-oxide layers 
of a (Insulator) / (LaMnOz)2 / (Sr MnOz)2 / (Insulator) het- 
erostructure for a — 4 A, e = 20, and C\ = 0.30. Figures are 
for (a) p(E ext = kV/cm); (b) Ap = p(E ext = 100 kV/cm) - 
p(E ext = kV/cm); and (c) M/pb at E ext = kV/cm and 
Eext = 100 kV/cm. MO layer 1 on the LMO side undergoes 
spin reversal when E ex t is applied. 



action energy, on the SMO side we take the e g electrons 
to be in one sub-lattice only (i.e., they are in the same 
sub-lattice in each layer but on alternate sub-lattices on 
adjacent layers). On account of virtual hopping, an e g 
electron polarizes all its neighboring sites that do not 
contain any e g electrons. Thus we observe that the half- 
filled interface will be fully polarized and that there will 
be a FM cluster that begins at the interface and perco- 
lates to the layers away from the interface on the SMO 
side. For instance, in the layer next to the oxide-oxide 
interface, all the e g electrons are next to sites unoccupied 
by e g electrons at the interface and hence are FM aligned 
with the interface; similarly, again in the adjacent layer, 
all the sites in the sub-lattice unoccupied by e g electrons 
have the same polarization as the sites occupied by e g 
electrons at the interface. The probability that a site, 
occupied (unoccupied) by an e g electron, contributes to 
the FM cluster is equal to the probability of finding the 
site occupied (unoccupied) multiplied by 1-P where P is 
the probability of not having any of the adjacent sites 
unoccupied (occupied) and belonging to the percolating 
cluster. Thus, we obtain the following set of coupled 
equations: 

a# = p(z N )[l-(l-2x u N _ 1 )(l-2x u N )\l-2x u N+1 )], (9) 



1 N 



0.5[1 - (1 - 2x^_ 1 )(l - 2.^) 4 (1 - 2a# , J], (10) 



where x£ O^jv) 1S the concentration of polarized sites that 
belong to the spanning cluster and these sites are a subset 
of the sites in the sub-lattice occupied (unoccupied) by 
e g electrons in the N th two-dimensional (2D) MO layer 
with TV being the index measured from the Insulator- 
SMO interface. Furthermore, z^ is the z-coordinate of 
the N th 2D MO layer. The factor (1 - 2x e £) [(1 - 2x u N )} 
represents the probability of a site that belongs to the 
sub-lattice occupied (unoccupied) by e g electrons in layer 
N but is not part of the spanning polarized cluster. The 
boundary conditions involving layer 1 are 

x e -=p{z l )[l-{l-2x^)\l~2x u 2 )], (11) 



xl = 0.5[1 - (1 - 2a; 1 9 ) 4 (l - 2x e 2 s )}, (12) 



Similar to the bulk situation^, in our heterostructure 
as well, we assume that in each MO layer two e g elec- 
trons on adjacent sites have the same ferromagnetic cou- 
pling J xy = 1.67 meV while two e g electrons on adjacent 
layers have again the same antiferromagnetic coupling 
J z = 1.2 meV. On the other hand, for an electron and a 
hole on adjacent sites in the same MO layer or adjacent 
MO layers, (due to virtual hopping of electron between 
the two sites) there is a strong ferromagnetic coupling 
J eh = t 2 /(2g 2 uJo) >> Jxy 1 ® ■ In our calculations, as long 
as the ratio of J xy /J z is fixed and J e h is the significantly 
dominant coupling, we get the same magnetic picture. 

When only a few MO layers are present in LMO, we 
will show that an LMO side that is not fully FM can be 
dramatically switched to being fully FM (by increasing 
the concentration of holes in the layers) upon application 
of a large external field. To establish this, we assume 
that the LMO-SMO interface and all MO layers up to 
M are completely polarized. Next, we assume that MO 
layers M and M — 1 have low density of holes so that 
there is a possibility that spins in MO layer M — 1 are 
not aligned with all MO layers up to layer M. We then 
analyze (in Appendix A) the polarization of MO layer 
M — 1 by comparing the energies for the following two 
cases: 

(i) layers M — 1 and M arc antifcrromagnctically 
aligned with the holes in layers M and M — 1 inducing 
polarization only on sites that are adjacent to the holes; 
and 

(ii) MO layer M — 1 is completely polarized and aligned 
with Layer M. 

In an LMO system that is three unit cell thick, for 
C\ = 0.93 eV (as shown in Fig. 2), we obtained a 
striking magneto-electric effect. For zero external field, 
when M = 2 is considered, case (i) has lower energy, 
i.e., layer 1 is antiferromagnetically coupled to layer 
2. On the other hand, when a strong electric field 
(~ 100 kV/cm) is applied, MO layer 1 (due to increased 
density of holes) becomes completely polarized and 
ferromagnetically aligned with the rest of the layers 
(i.e., MO layers 2, 3, and the oxide-oxide interface). 
Thus, we get a giant magneto- electric effect! Similarly, 
in a {Insulator) / {LaMn03)2 / (SrMnOz)2 / (Insulator) 
heterostructure, for C\ = 0.30 eV (as shown in Fig. 3), 
we again obtain a huge magneto-electric effect. 



while those for the LMO-SMO interface are 
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CONCLUSIONS 



B. LMO side 

Next, we will show that a few-layered LMO side can 
be border-line ferromagnetic [i.e., it is either barely fully 
ferromagnetic or close to being fully ferromagnetic] which 
can be exploited to obtain giant magneto-electric effect. 



We used the fact that when the two Mott insulators 
LMO and SMO are brought together to form a het- 
erostructure, one may realize the entire phase diagram 
of La\- x Sr x MnO^ across the heterostructure with the 
x = phase occurring at the LMO surface at one end and 
evolving to the x — 0.5 state at the oxide-oxide interface 
and finally to the x = 1 phase at the other end of SMO. 



By changing the charge profile using an applied electric 
field, we showed that the FMI region can be expanded at 
the expense of the A-type AFM domain thereby produc- 
ing a giant magnetoelectric effect. Since the competing 
phases at the interface of the FMI and the A-type AFM 
region are close in energy, the interface can be tuned to 
be on the verge of a FM transition by the application 
of a magnetic field of the order of a Tesla (as the Zee- 
man energy is comparable to the difference in the FM 
coupling J xy and the AFM coupling J z )\ then, applying 
a tiny electric field will produce a FM transition at the 
interface. 

We also would like to mention that if a su- 
perlattice were formed from the heterostructure 
(Insulator) / \LaMnOs) n / (SrMnOs) n / (Insulator) 1 
then the dipoles from each repeating heterostructure 
unit will add up to produce a giant electric dipolc 
moment. On the other hand, if the insulating layers 
were not present, the superlattice formed from the 
repeating unit (LaMnOz) n /(SrMnO^) m would not 
produce a large electric dipole as the charge from LMO 
can leak to SMO on both sides leading to a small net 
dipolc moment. 

Lastly, minimizing the energy expressed as a functional 
of density and magnetization at a fixed electric field (so 
as to get simultaneously both charge and magnetization 
profiles) is more accurate (albeit tedious) than solving, 
as we did, for the charge profile and using it to get the 
magnetization picture. This better alternative can also 
address the issue of how an external magnetic field affects 
the density profile and the ferroelectric moment. 



Appendix A 

In this appendix, we will analyze the nature of mag- 
netism on the LMO side. To begin, we assume that all 
MO layers from M to the oxide-oxide interface (i.e., MO 
layers M, M + 1, ...) arc completely polarized. We then 
bring a completely polarized MO layer to form the layer 
M — I next to the MO layer M and wish to determine 
the interaction energy between layer M — 1 and layer M 
and their relative polarization. We compare the energies 
in the following two scenarios: 

(i) MO layer M — 1 is antiferromagnetically aligned with 
layer M and there is an induced polarization next to the 
holes in these layers, i.e., all the sites that are adjacent 
to the holes in layers M and M — 1 end up having the 
same polarization as the holes; and 
(ii) MO layer M — 1 is completely polarized and ferro- 
magnetically aligned with Layer M . 

The energy of interaction in case (i) is the sum of the 
two terms T a and XJ, obtained below: 
(a) when an electron in the electron sub-lattice of 



layer M is not surrounded by a hole in either layer 
M or layer M + I [the probability of which is 

(I — 2p(zm)) (1 — 2p(zM+i))], then the hole in layer 
M — I flips the spin of the electron below it in layer 
AI. Contrastingly, a hole in layer M will always flip the 
spin of the electron above it in layer M — 1. The resul- 
tant electron- hole interaction energy is —J e h x [p(zm) + 
p(zm-i)\- The flipping of spins by holes changes the in- 
teraction between the spins and their four neighbors in 
the xy-plane and yields the interaction energy 2 x J xy x 
[p(zm-i) x 4 + p(zm) x 4]. Furthermore, there is an ad- 
ditional antiferromagnetic coupling between layer AI and 
AI — 1 which is given by — J z x [I — p(zm)] where the anti- 
ferromagnetic coupling between layers is reduced due to 
holes in layer M . Thus the resultant interaction energy 
is 



T a = 



(1 - 2p(z M )) (1 - 2p(zm + i)) 

X [-Jeh X \p(z M ) + p(ZM-l)} 

+ 2J xy x [Ap(z M -i) + 4p(z M )} 
-J z x [1 - p(z M )}} ■ 



(Al) 



(b) When an electron in the electron sub-lattice of 
layer M is surrounded by at least one hole in either 
layer M or layer M + 1 [the probability of which is 
1 — (I — 2p(zu)) (1 — 2p(zA/+i))], then the hole in layer 
AI — I and its 4 neighbors in the same layer flip their 
spins to align with the spin of the electron below the 
hole in layer AI — 1. On the other hand, a hole in layer 
M always flips the spin of the electron above it in layer 
AI — I. The resultant electron-hole interaction energy is 
— J e h x [p(zm) + p(zm-i)]- The flipping of spins changes 
the interaction between the flipped spins and their neigh- 
bors in the plane and yields the interaction energy 
2xJ xy x [p(z M ~i) x4x3 + p(zm)x4]. Additionally, there 
is a new antiferromagnetic coupling between layer M and 
AI -1 given by -J z x [1 - p(z M ) - p(z M -i) x (1 + 4x2)] 
where the antiferromagnetic coupling between layers is 
reduced due to holes in layers M and AI — 1 and flipped 
electron spins in layer AI — I. Thus the resultant inter- 
action energy is 



T h 



1 - (1 - 2p(z M )) 4 (1 - 2p(z M+1 )) 

' [-Jeh X [p{z M ) + p(z M -l)} 

+ 2J xy x [\2p(z M -i) + M z m)) 
-J z x [1 - p(z M ) - 9p(z M -i)}} ■ 



(A2) 



Next, the interaction energy for case (ii) is given by 
sum of the following two terms T c and X^: 

(c) T c = -Jeh x [p(z M ) + p(z M -i)}; and 

(d) T d = J z x [1 - p(z M ) - p(z M -i)}- 
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